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Abstract. The Cold Dark Matter paradigm successfully explains many 
phenomena on scales larger than galaxies, but seems to predict galaxy ha- 
los which are more centrally concentrated and have a lumpier substructure 
than observed. Endowing cosmic dark matter with a small primordial veloc- 
ity dispersion preserves the successful predictions of the Cold Dark Matter 
scenario on large scales and improves the agreement with halo structure. 
A "phase density" Q, proportional to the inverse entropy for nonrelativis- 
tic matter, is estimated for relativistically decoupled thermal or degenerate 
relic particles of mass mx, with a numerical factor depending on the par- 
ticle type but no cosmological parameters. Since Q cannot increase for 
dissipationless, collisionless matter, at a given velocity dispersion there is 
a maximum space density; this "phase packing" constraint eliminates the 
singular density predicted by CDM. The core radius and halo circular ve- 
locity scale analogously to degenerate dwarf stars. Particle velocities also 
filter primordial perturbations on a scale depending on Q and on details of 
particle distributions. Particle candidates for warm matter are briefly dis- 
cussed; for warm thermal relics to have the observed mass density requires 
decoupling prior to the QCD epoch and therefore a superweak interaction 
with thermal Standard Model particles. 
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1. How Cool is Cold Dark Matter? 

The successful match of predictions for large scale structure and microwave 
anisotropy vindicates many assumptions of standard cosmology, in particu- 
lar the hypothesis that the dark matter is composed of primordial particles 
which are cold and collisionless. At the same time, the CDM paradigm 
finds difficulty explaining the small-scale structure within galaxy haloes: 
CDM appears to predict excessive relic substructure [1, 2] in the form of 
many dwarf galaxies which are not seen and may disrupt disks, and also 
predicts[3, 4, 5, 6, 7] a universal, monotonic increase of density towards 
the center of all halos which is not seen in close studies of dark-matter- 
dominated galaxies. The latter problem seems to arise as a generic result of 
low-entropy material sinking during halo formation, quite independently of 
initial conditions; the former effect seems to be a generic result of hierarchi- 
cal clustering predicted by CDM power spectra, which produce fluctuations 
on small scales that collapse early and survive as substructure. Although 
these problems are still controversial from both a theoretical and obser- 
vational point of view [8], it is not easy to dismiss these effects by various 
complicated baryonic devices. 

It is possible that the problems with halo structure are giving specific 
quantitative clues about new properties of the dark matter particles. The 
existence of dwarf cores and smooth substructure are just what one would 
expect if the dark matter is not absolutely cold but has a small nonzero 
primordial velocity dispersion. Such a model produces two separate effects: 
a phase packing or Liouville limit which produces halo cores, and a filter 
in the primordial power spectrum which limits small-scale substructure. 
The estimated dispersion required for the two effects does not quite agree 
but is close enough to motivate a closer look. This discussion is meant to 
motivate more detailed comparison of models and theory aimed at using 
halo properties to test the hypothesis of primordial velocity dispersion and 
ultimately to measure particle properties. 

The physics of the filtering by freely streaming particles closely parallels 
that of massive neutrinos[9, 10, 11, 12, 13], the standard form of "hot" dark 
matter. A thermal particle which is more weakly interacting and therefore 
separates out of equilibrium earlier than neutrinos, when there are more 
particle degrees of freedom, has a lower temperature than neutrinos and 
therefore both a lower density and lower velocity dispersion for a given 
mass. Such "thermal relics" constitute one class of warm dark matter can- 
didate; there are other possibilities, including degenerate particles[15, 16] 
and products of decaying particles. In many respects their astrophysical ef- 
fects are very similar to the thermal case since up to numerical factors the 
damping scale and phase packing limit are both fixed by the same quantity, 
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the classical "phase density" of the particles. For particles which separate 
out when relativistic, the phase density depends only on the particle prop- 
erties and not on any cosmological parameters. It is similar (though not 
identical) for bosons and fermions, and in thermal and degenerate limits. 

The phase-packing constraint is also familiar from the context of massive 
neutrinos. Tremaine and Gunn[17] showed that the phase density of dark 
matter in giant galaxies implies a large neutrino mass, therefore too large 
a mean cosmic density. However, their argument can be turned around to 
explain the lack of a cusp at the center of dwarf halos. An upper limit to 
the central density of an isothermal sphere can be derived for a given phase 
density; a very rough comparison with dwarf dynamics suggests a limit 
corresponding to that of a thermal relic with a mass of about 200 eV. This 
is lighter (that is, a larger dispersion) than the «1 keV currently guessed 
at from the filtering effect, but not by so much that the idea should be 
abandoned; the simulations and comparison with data may yet reconcile 
the two effects and could reveal correlations of core radius and velocity 
dispersion predicted by the phase-packing hypothesis. Comparison of the 
two effects may also reveal new dark matter physics: for example, if the 
particles scatter off each other by self-interactions (which may even be 
negligible today), free streaming is suppressed and the filtering occurs on 
a somewhat smaller scale. 

Warm dark matter has most often been invoked as a solution to fixing 
apparent (and no longer problematic[18]) difficulties with predictions of the 
CDM power spectrum for matching galaxy clustering data. A filtered spec- 
trum may however solve several other classic problems of CDM on smaller 
scales, in galaxy formation itself. Baryons tend to cool and collapse early 
into lumps smaller and denser than observed galaxies [19]. Although this 
may be prevented by stellar-feedback effects, recent studies[20] suggest that 
CDM has fundamental difficulties explaining the basic properties of galax- 
ies such as the Tully-Fisher relation; dynamical loss of angular momentum 
results in halos which are too concentrated. Some of these problems may be 
solved[21] in warm dark matter models which suppress the early collapse 
of subgalactic structures. Modeling disk formation includes baryonic evolu- 
tion so requires understanding the ionization history of gas; an important 
constraint comes from the observed structures of the Lyman-a forest [22]. 
Simulations[21] suggest that the optimal filtering scale corresponds to a 
thermal particle mass of about 1 keV. 

For these relatively massive thermal relics to have the right mean den- 
sity today, the particle must have separated out at least as early as the 
QCD era, when the number of degrees of freedom was significantly larger 
than at classical weak decoupling. Its interactions with normal Standard 
Model particles must be "weaker than weak," ruling out not only neutrinos 
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but many other particle candidates. The leading CDM particle candidates, 
such as WIMPs and axions, form in standard scenarios with extremely high 
phase densities, although more elaborate mechanisms are possible to endow 
these particles with the required velocities. It is therefore of considerable in- 
terest from a particle physics point of view to find evidence for the existence 
of finite primordial phase density from galaxy halo dynamics. Neither the 
theory or the observational side allow a definitive case to be made as yet, 
but the evidence is certainly suggestive. In principle, dynamics can provide 
detailed clues to the dark matter mass and interactions. Here I give a few 
examples of simple calculations which reveal the connections between the 
particle properties and the halo properties. 



2. Phase Density of Relativistically-Decoupled Relics 

Consider particles of mass m originating in equilibrium and decoupling 
at a temperature To » m or chemical potential \i » m. The original 
distribution function is [23] 

/(p) = ( e (*H0/T D ± 1)-1 ~ (e (p-M)/T D ± (1) 

with E 2 = p 2 + m 2 and ± applies to fermions and bosons respectively. The 
density and pressure of the particles are [24] 

where g is the number of spin degrees of freedom. Unless stated otherwise, 
we adopt units with h = c = 1. 

With adiabatic expansion this distribution is preserved with momenta 
of particles vary as p oc so the density and pressure can be calculated 
at any subsequent time[25]. For thermal relics fx = 0, we can derive the 
density and pressure in the limit when the particles have cooled to be 
nonrelativistic: 

n=^-(^- (4) 
(2vr)3 J eP ± 1 1 j 

5 T 5 f p 2 d 3 p ^ 



(2vr) 3 3m7 e?> ± 1 

where the pseudo-temperature To = Td(Rd/Ro) records the expansion of 
any fluid element. 
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It is useful to define a "phase density" Q = p/(v 2 ) 3 ^ 2 proportional 
to the inverse specific entropy for nonrelativistic matter. This quantity is 
preserved under adiabatic expansion but for nondissipative particles cannot 
increase. Combining the above expressions for density and pressure and 
using (v 2 ) = 3P/nm, we find 

Qx = qxgxm x . (6) 



The coefficient for the thermal case is 

y (2vr) 3 [Jdp(p i /eP± l)]3/2 ' v ; 

where the last equality holds for thermal fermions. An analogous calcula- 
tion for the degenerate fermion case (T = 0, [Id >> fnx) yields the same 
expression for Q but with a different coefficient, 

q^-^ 11 */^ = 0.036335. (8) 

(27T)3 [ /o V dp] 3/2 ^ 

The phase density depends on the particle properties but not at all on 
the cosmology; the decoupling temperature, the current temperature and 
density do not matter. Up to numerical factors (which depend on thermal 
or degenerate, boson or fermion cases), the phase density for relativistically 
decoupled or degenerate matter is set just by the particle mass mj. 



3. Phase Packing and the Core Radius of an Isothermal Halo 

For a given velocity dispersion at any point in space, the primordial phase 
density of particles imposes an upper limit on their density p. Thus dark 
matter halos do not form the singular central cusps predicted by Cold Dark 
Matter but instead form cores with constant density at small radius. A lower 
limit to the size of the core can be estimated if we assume that the matter 
in the central parts of the halo lies close to the primordial adiabat defined 
by Q. This will be good model for cores which form quietly without too 
much dynamical heating. This seems to be the case in typical CDM halos, 
indicated by the cusp prediction; it could be that warm matter typically 
experiences more additional dynamical heating than cold matter, in which 
case the core could be larger. 

The conventional definition of core size in an isothermal sphere [26] is 
the "King radius" 

r = ^9a 2 /47rGp (9) 
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where a denotes the one-dimensional velocity dispersion and p denotes the 
central density. Making the adiabatic assumption, po = <2(3cr 2 ) 3 / 2 , we find 

r = yWi/^3^2(QG Wcoo )- 1 / 2 = 0.44(QG^ COO )- 1 / 2 (10) 

(Note that aside from numerical factors this is the same as a degenerate 
dwarf star; the galaxy core is bigger than a Chandrasekhar dwarf of the 
same specific binding energy by a factor (m pro ton/ f^x) 2 ■) For comparison 
with observations we have expressed the core radius as a function of the 
asymptotic circular velocity v coo = y/2a of the halo's flat rotation curve. 
For the thermal and degenerate phase densities derived above, 

r ,thermal = 5.5kpc(m x /100eV)- 2 (t; coo /30kms- 1 )- 1 / 2 (11) 

^degenerate = 1.3kpc(mx/100eV) " 2 (WSOkmS^ 1 )~ V2 , (12) 

where we have set g = 2. The circular velocity in the central core displays 
the harmonic behavior v c oc r; it reaches half of its asymptotic value at a 
radius r\/2 ~ 0.4ro. 

The best venue for studying the effect is in the small, dark-matter- 
dominated disk galaxies where the cusp problem of CDM seems most clearly 
defined. Material on circular orbits gives a direct measure of the enclosed 
mass and therefore of the density profile. A rough guess from current ob- 
servations of inner rotation curves of a few dwarf spiral galaxies suggests a 
halo core corresponding to a thermal particle mass of about 200 eV[27]; for 
a larger mass, additional (nonprimordial) dynamical heating is required. 
This is also consistent with what is known about dark matter in dwarf 
elliptical galaxies from studies of stellar velocities[28, 29, 30]. 

The relationship of core radius with halo velocity dispersion is a simple 
prediction of the primordial phase density explanation of cores which will 
probably generalize in some form to a cosmic population of halos. In par- 
ticular if phase packing is the explanation of dwarf galaxy cores, the dark 
matter cores of giant galaxies and galaxy clusters are predicted to be much 
smaller than for dwarfs, unobservably hidden in a central region dominated 
by baryons. 

4. Filtering of Small-Scale Fluctuations 

The transfer function of Warm Dark Matter is almost the same as Cold 
Dark Matter on large scales, but is filtered by free-streaming on small 
scales. The characteristic wavenumber for filtering at any time is given by 
kx ~ H/(v 2 ) 1 / 2 , with a filter shape depending on the detailed form of the 
distribution function. In the current application, we are concerned with H 
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during the radiation-dominated era [z > 10 4 ), so that H 2 = 8TtGp re i/3 oc 
(l + z) 2 , where p re i includes all relativistic degrees of freedom. For constant 
Q, (f 2 ) 1 / 2 = (px/Q) 1 ^ 3 (1 + z ) as long as the X particle are nonrela- 
tivistic. The maximum comoving filtering scale[24] is thus approximately 
independent of redshift and is given simply by 

kx « H^v- X \ (13) 

where £t re i = 4.3 x 10~ 5 /i~ 2 and vxo = (Q/pxo)^ 1 ^ 3 is the rms velocity of 
the particles at their present mean cosmic density pxo'- 

k x ~ 0.65Mpc _1 (u X o/lkm s" 1 ) -1 , (14) 

with no dependence on Hq. For the thermal and degenerate cases, in terms 
of particle mass we have 

vxo,thermal = 0.93km s" 1 h% 3 \m x /lOOeV)" 4 / 3 (Q x /0.3) 1 / 3 <T 1/3 (15) 

V X0 ,degenerate = 0.35km S - 1 /l^ 3 (m X /100eV)- 4 / 3 (^ X /0.3) 1 / 3 5 - 1 / 3 (16) 

In the case of free-streaming, relativistically-decoupled thermal parti- 
cles, the transfer function has been computed precisely [12, 21]; the charac- 
teristic wavenumber where the square of the transfer function falls to half 
the CDM value is &a/2, stream = k x /5.5. The mass implied for this kind of 
candidate to preserve the success of CDM on galaxy scales and above is 
about[21] 1 keV; if it is much smaller (in particular, as small as the 200 
eV we require for phase packing alone to help the core problem), filtering 
occurs on too large a scale. (A filtering scale of roughly k pa 3/i7oMpc _1 
preserves the successes of CDM on large scales and helps to solve the CDM 
predictions of excess dwarf galaxies, excessive substructure in halos, in- 
sufficient angular momentum, and excessive baryon concentration.) This 
problem might be fixed in other models with a different relationship of £4/2 
and kx- For example, if the particles are self-interacting the free streaming 
is suppressed and the relevant scale is the standard Jeans scale for acoustic 
oscillations, kj = y/3H/cs = \/27/5kx, which is significantly shorter than 
k\/2,stream a * a fixed phase density. Alternatively it is possible that warm 
models might be more effective at producing smooth cores than we have 
guessed from the minimal phase-packing constraint above; an evaluation of 
this possibility requires simulations which include not just a filtered spec- 
trum but a reasonably complete sampling of a warm distribution function 
in the particle velocities[31]. 

5. Density of Thermal Relics 

A simple candidate for warm dark matter is a standard thermal relic — a 
particle that decouples from the thermal background very early, while it 
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is still relativistic. In this case the mean density of the particles can be 
estimated [24] from the number of particle degrees of freedom at the epoch 
T D of decoupling, g* D : 

n x = 7.83h~ 2 lg eff /g* D ](m/W0eV) = 0.24/ l70 2 (m x /100eV)( 5e// /1.5)(^ Z) /100)- 1 

(17) 

where g e f f is the number of effective photon degrees of freedom of the 
particle (= 1.5 for a two-component fermion). For standard neutrinos which 
decouple at around lMeV, g*D = 10.75. 

An acceptable mass density for a warm relic with rax > 200 eV clearly 
requires a much larger g*D than the standard value for neutrino decoupling. 
Above about 200 MeV, the activation of the extra gluon and quark degrees 
of freedom (24 and 15.75 respectively including uds quarks) give g*o ~ 50; 
activation of heavier modes of the Standard Model above ~ 200GeV pro- 
duces g^d ~ 100 which give a better match for fix < 0.5, as suggested 
by current evidence. Masses of the order of 1 keV can be accomodated by 
adding extra, super symmetric degrees of freedom. Alternatively a degener- 
ate particle can be introduced via mixing of a sterile neutrino, combined 
with a primordial chemical potential adjusted to give the right density [16]. 
Either way, the particle must interact with Standard Model particles much 
more weakly than normal weak interactions, which decouple at ~ 1 MeV. 

Note that all warm dark matter particles have low densities compared 
with photons and other species at 1 MeV so they do not strongly affect 
nucleosynthesis. However, their effect is not entirely negligible. They add 
the equivalent of (Tx/T^) 3 = 10.75/ g*d ~ 0.1 to 0.2 of an effective extra 
neutrino species, which leads to a small increase in the predicted primordial 
helium abundance for a given rj. Because the phase density fixes the mean 
density at which the particles become relativistic, this is a generic feature for 
any warm particle. This effect might become detectable with increasingly 
precise measurements of cosmic abundances. 
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